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Abstract 

The kinematics of a final state system with two invisible particles and two visible particles can 
develop cusped peak structures. This happens when the system has a fixed invariant mass (such 
as from a narrow resonant particle decay or with a fixed collision cm. energy) and undergoes 
decays of two on-shell intermediate particles. Focusing on the "antler decay topology", we derive 
general analytic expressions for the invariant mass distribution and the kinematic cusp position. 
The sharp cusp peaks and the endpoint positions can help to determine the masses of the missing 
particles and the intermediate particles. We also consider transverse momentum variables and 
angular variables. In various distributions the kinematic cusp peaks are present and pronounced. 
We also study the effects on such kinematic cusp structures from realistic considerations including 
finite decay widths, the longitudinal boost of the system, and spin correlations. 
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I. INTRODUCTION 



With the advent of the Large Hadron Collider (LHC), the TeV scale physics will be 
fully explored in the coming decades. Most pressing of all to learn is the mechanism of the 
electroweak symmetry breaking and the related underlying dynamics beyond the standard 
model (SM). Among many interesting phenomena associated with the new physics at the 
TeV scale, the signature of events with large missing energy is one of the most exciting 
possibilities at the LHC. This is expected from new particles that do not leave any trace in 
the hadronic and electromagnetic components of the detector. These new missing particles 
may help to address one of the most profound puzzles in cosmology: what constitutes 
nearly a quarter of the energy density of our current universe in a form of cold dark matter 
(CDM) The thermal history of the early universe suggests that a stable neutral particle of 
the electroweak-scale mass and interaction, called the Weakly Interacting Massive Particles 
(WIMP), is a plausible explanation of CDM [2] and may be discovered as a missing particle 
at TeV-scale colliders. 

Missing energy signal is generic in many new physics models. Additional discrete sym- 
metry is often introduced to prohibit dangerous processes such as proton decay and to make 
the model compatible with the electroweak precision tests. Such a discrete symmetry (or 
parity) often needs nontrivial representations of new particles, while it assigns vanishing 
charges (or trivial representation) to the SM particles. Therefore, the lightest new particle 
is stable, becoming a natural candidate for the CDM particle. One of the most studied 
examples is the lightest neutralino in supersymmetric (SUSY) theories with i?-parity con- 
servation [3]. Other examples include the lightest Kaluza-Klein (KK) particle in universal 

extra dimensional (UED) theories with KK parity conservation (4], and the heavy photon in 

r~i 

the little Higgs models with T-parity conservation (LHT) [5j. In this regard, the search for 
missing particles at the LHC and future colliders has great implications in understanding 
both the fundamental particle physics and the nature of our universe. At hadron colliders, 
the experimentally observable signature will be missing energy-momentum transverse to the 



beam direction. Great efforts have been made on t 
ing energy signals in various new physics models 



re phenomenological studies of the miss- 
7J and optimistic conclusions have been 



reached such that significant excess is expected above the SM background at the LHC js]. 
In order to reveal the CDM identity and to compare with the results from direct and 
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indirect DM searches, it is imperative to determine the missing particle mass at colliders. 
This is a very challenging task since such weakly interacting neutral particles leave neither 
charged particle tracks nor significant energy deposit in the detector. Furthermore, the 
missing particles always come in pairs in an event due to the conserved "parity" , so that the 
final state kinematics is under-constrained. Finally, if we consider hadronic collisions such as 
at the LHC, the partonic cm. energy as well as the frame are unknown on an even-by-event 
basis. 

As reviewed in Ref. Jq], most of the techniques for the missing particle mass measurement 
at the LHC can be categorized into the following three cases: (i) endpoint methods 101; 

nfi nn 

(ii) polynomial methods [UJ, ll2j; (iii) M^ 2 methods |13Ml6l|. All three methods rely on a 



cascade decay of a heavy new particle, ended up with a single missing particle X. At each 
step of a cascade chain, a visible particle is produced, which may provide information on 
the missing particle mass as well as the intermediate new particle mass. 

Endpoint methods use the kinematic edges of invariant mass distributions of the visible 
particles in a given cascade decay. If the cascade chain is long enough with at least three 
visible particles, the number of kinematical constraints is sufficient to determine all the 
masses involved [9|. When the decay chain is not long, the observables are insufficient for 
complete mass determination. In addition, the positions of endpoints are more sensitive to 
the mass difference than to the absolute mass. 

Polynomial methods use reconstructable events in which the number of the on-shell kine- 
matic constraints exceeds the number of the unknown masses and momentum components. 
By combining multiple event information, one can maximize the information for determi- 

n 

nation of mass parameters However, this method intrinsically requires a long decay 

chain, at least two-step cascade decays in each chain, producing four visible particles {9]. It 
suffers from small statistics and large combinatoric background. 

The Mt2 variable, originally proposed in Ref. jl3|, is useful at hadron colliders for mea- 
suring the mass of a new mother particle when pair-produced. Two mother particles decay 
through the same decay chain. For each chain, the transverse mass is constructed with the 
missing transverse momentum. As a function of a trial mass for the missing particle, My 2 is 
the minimum value of the larger value of these two transverse masses. The minimization is 
over all possible missing transverse momenta of two decay chains as satisfying the observed 
total missing energy constraint. The M?2 distribution has the maximum at the mother par- 
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FIG. 1: The antler decay topology of a heavy new particle D into two missing particles (X\ and 
X2) and two visible particles {a\ and 02). 
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FIG. 2: The cascade decay topology of a heavy new particle D into two missing particles {X\ and 
X2) and two visible particles {a\ and 02). 

ticle mass when the trial mass hits the true missing particle mass. Therefore, it provides one 
relation between the mother particle mass and the missing particle mass. A more exciting 
observation is that the endpoint curve of Mrp, as a function of the trial mass shows a kink 



where the trial mass becomes the true mass 



1$. 



In all three methods above, a crucial issue is how to fully reconstruct the kinematics of a 
signal event. This relies on exclusive selection of events of a given type. If the decay chain is 
long, the reconstruction becomes more difficult as combinatoric complications emerge: the 
large number of involved particles entangle the origin of the decay of each observed particle. 
The hemisphere method, an algorithm to group collinear and high-p^ particles, was shown 
to be useful to some extent in the inclusive Mt2 analysis for the disentanglement of the 
data Q- 

Recently, it has been pointed out that the missing particle mass can be determined from 



singular structures in kinematic distributions for shorter and simpler decay chains 



18 



19] 



In our previous work [18J, we considered a resonant "antler decay" of a heavy new particle 
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into a pair of missing particles and a pair of SM visible particles, as shown in Fig. [TJ and 
found non-smooth peaks in some kinematic distributions. These peaks are called "cusps" 
and the positions of the cusps depend only on the masses of the involved particles. The cusp 
is statistically more advantageous because it is at the peak region. The mass measurements 
can be benefited from knowing the kinematic cusp structures. 

We consider a resonant decay of a heavy particle D into two visible particles and two 
missing particles. The invariant mass distributions of this type of decay were first presented 



in our previous publication |18|, and recently further studied in Ref. j^j. Obviously this 



heavy particle D is parity-even. The general topology of such resonance decays is divided 
into two classes: 

1. Antler decay topology: a heavy particle D decays into two parity-odd particles {B\ and 
B 2 ) at the first step and each parity-odd particle subsequently decays into a missing 
(denoted by dashed lines) particle and a visible particle, as in Fig. [TJ 

2. Cascade decay topology: a heavy particle D splits to two particles with one or both 
visible particles at each step, finally into a missing particle. According to at which 
step the first missing particle comes out, there are two non-trivial cascade topologies, 
as in Fig. [2j 

The antler decay and the cascade decay are siblings to each other as they share the same 
skeleton of topology. Since they have different orientation of incoming and outgoing particles, 
the cusps appear with different manifestations. In this paper, we focus on the antler top olog y 



2l|. 



only and leave the presentation on the cascade decay topology to a companion paper 

Antler decays arise in many new physics models. We now list a few examples for illus- 
tration. 

re heavy CP-even neutral 



In the Minimal Supersymmetric Standard Model (MSSM), t 
Higgs bosons may have sizable rates of the following decay 221 ] : 

H^xl + xl^ZxUZXv (1) 

In the MSSM with an additional U(l) gauge interaction, the extra U(l) gauge boson 



Z' can have antler decay modes like 23] 



z'^r + e+->e-xl + e + xv (2) 
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The ultraviolet completion of the LHT model often involves an extension of the Higgs 
sector that accommodates heavy Higgs bosons. Large top Yukawa coupling leads to 
substantial decay of the neutral heavy Higgs into a pair of T-parity odd top quarks 
t_, followed by t_ decay into the SM top quark and the heavy photon Ah (the CDM 
candidate) 

H ->■ t_ + t_ ->■ tA H + iA H . (3) 
In the UED model with KK parity conservation, the second KK mode of the Z boson 



can have antler decay modes 24j. Z^ decays into a pair of the first KK modes of the 



lepton, followed by its decay into a SM lepton and the CDM particle 

Z V> -> L« + tfiW + tBW. (4) 

• At lepton colliders with e + e~ or collisions, the well-determined cm. energy 

renders some pair production and their subsequent decay processes to be of the antler 
topology. One example is 

e + e-/fi + fi- -»• £ + + 1 -> £ + Xi + £~Xi- (5) 

In the current work, we only focus on the generic features of antler kinematics. The rest of 
the paper is organized as follows. We begin our discussion by explaining the unique features 
of the antler kinematics in Sec. [Til Focused on the symmetric antler decay, we consider the 
massive visible particle case, and demonstrate in Sec. II III the cusps and endpoints in the 
kinematic distributions of the invariant mass, transverse momenta, and angular variables, 
constructed from two visible particles. In Sec. IIV[ we study the massless visible particle 
case. We discuss some effects of more realistic considerations in Sec. El such as the finite 
decay widths of the resonant particles, the longitudinal boost between the cm. frame and 
the lab frame, and spin correlations. We conclude in Sec. I VI I A few appendices are devoted 
to some technical details for a general four-body phase space treatment, the derivations of 
the cusp peak and analytic expressions of some kinematic distribution, and more discussions 
for the general antler decay. 
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FIG. 3: Kinematic configuration of the antler decay in the rest frame of the parent particle D via 
two intermediate particles B\ and B2, followed by B — > aX. 9\ and 82 are defined in the rest 
frames of B\ and B2, respectively, and <j) in the D rest frame. 

II. KINEMATICS OF ANTLER DECAY TOPOLOGY WITH TWO MISSING 
PARTICLES 

We consider the resonant decay of a heavy particle D into two visible particle a± and 02, 
and two missing particle X\ and X2 via a chain of two-body decays through intermediate 
particles Pi and B 2 , as depicted in Fig. [TJ 

D{P) P 1 (p 1 )+P 2 (p 2 ), (6) 
Pi(pi) -> ai(A;i) + Xi(fc 3 ), P 2 (p 2 ) -»■ 02(^2) + X 2 (k 4 ). 

Since most of the processes of our interest are symmetric between two decay branches, we 
focus on the symmetric antler decay, defined by 

Symmetric antler: B = By — B 2 , a = a\ = a 2 , X = X± = X 2 . (7) 

The general antler decay with arbitrary masses is to be discussed in Appendix [Bj 

In the three-dimensional momentum space, the kinematic configuration of the antler 
decay is illustrated in Fig. [3j In the rest frame of the parent particle D, the intermediate 
particles B\ and B 2 are moving back-to-back, and the momentum direction defines the 
principal decay axis z, with B\ moving into the +z direction and B2 into the —z direction. 
Two momenta of a\ and X\ in the B\ rest frame form the decay plane PI, which is identified 
as the xz-plane. In the same way, the decay plane PI is defined by the B 2 decay products. 
In the decay plane PI, we define a polar angle 61 between the +z direction and the a\ 
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momentum in the B\ rest frame. Similarly, 62 is the polar angle between the —z direction 
and the a 2 momentum in the B2 rest frame. The azimuthal angle between two decay planes 
PI and P2 is denoted by 0. 

As explicitly shown in Appendix these three internal angles (61, 8 2 , and <fr) represent 
the phase space configuration of the antler decay topology. The dynamics of the antler decay 
is encoded in the differential decay width dT defined in the rest frame of D. dT is a function 
of the internal phase space variable (61, 62, <$)'■ 

dT oc \M\~ 2 d$ 4 , ( 8 ) 

where \Ai\ 2 is a reduced matrix elements and <i$4 = d cos 61 dcos6 2 d<p (see Appendix lAl for 
more details). The reduced matrix element \A4\ 2 is a smooth function of (61, 82, <p), and thus 
dT/d^4 does not show any singular behavior. 

Kinematic singularities emerge as missing particles allow us only the projection of the 
full kinematic phase space onto a lower dimensional phase space accessible by the visible 
particle momenta. This partial access inevitably hides some of necessary information for 
the full mass reconstruction. However, we can still decode the mass information out of some 
observables, say F's. 

In order to obtain dT/dY, we project the hypersurface of the phase space (9i,9 2 ,(f)) 
onto Y: for each value of Y, dT/dY is proportional to the volume of the hypersurface 
corresponding to that specific value of Y. When the hypersurface fails to be a manifold at a 
certain point Y, dT/dY develops non-smoothness. This is called singular points, where the 
differential (dY/d8i, dY/d8 2 , dY/dcf)) vanishes 1 . As a result, we see non-smooth behaviors in 
the distribution of Y, which give rise to kinematic cusps and endpoints. General discussions 
on the development of singularity in the multi-dimensional observable phase space have been 



presented in Ref. |l9] 



Since the parent particle D is moving in the lab frame, the observable variable Y from 
the momenta of two visible particles a\ and a 2 can be classified into three categories: 

• Lorentz-invariant observable: there is only one Lorentz-invariant observable, the in- 



1 In multi-dimensional cases, this condition is a reduced rank condition of Jacobian matrix of mapping from 
the phase space to the observable Y's 19 1. 
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variant mass of a\ and a 2 , 

m = V(h + k 2 ) 2 . (9) 

Longitudinal-boost invariant observables: 

— the transverse momentum of one visible particle i: 

PTi = |k; T | . (10) 

Here and henceforth, a bold-faced letter denotes a three-momentum. 

— the total transverse momentum of the a\-a 2 system: 

p T = |k 1T + k 2T \ ■ (11) 

In the four-body decay under consideration, this is the same as the magnitude of 
the missing transverse momentum j>T of the decay. 

— the transverse mass of the a\-a 2 system: 



rriT = y Pt + m2 ■ (12) 

— the cluster transverse mass of the a\-a2-j>T system: 

m% = rax + fir- (13) 

— the rapidity difference: 

A?7 = \rj a i -r] a2 \, where r] ai = ^ In ^ k ^ z ^j . (14) 

Non-invariant observable: we consider an angular variable 0, which is the angle be- 
tween one visible particle (say ai) in the cm. frame of a\ and a 2 and the cm. moving 
direction in the D rest frame, given by 

kj" ia2) -(k 1 + k 2 )^> 
cos6 = -— — -Tjrc. (15) 

|kS aia2) | |k 1 + k 2 | (D) 

In what follows, the superscript in a momentum denotes the reference frame. In the 
main text, a momentum without a superscript is in the lab frame. 
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III. MASSIVE VISIBLE PARTICLE CASE 



In this section, we consider the case of massive visible particles. For a resonant decay, it 
is very convenient to express the kinematics in terms of rapidity variables. For a two-body 
decay of % — \ j + k, we write the four-momentum of the particle j in the rest frame of the 
parent particle i as Pj = (Ej , p^) = (rrij coshr^-, rrij sinhrjjPj). Here rjj is the rapidity of 
particle j in the rest frame of the parent i, given by 



E 

coshr/j = 



(i) 



rrij 



mf + rrij — ml 
2rriimj 



(16) 



The superscript of a rapidity, specifying the reference frame, is omitted when it is the rest 
frame of the parent particle. In this section, we assume that all the particles are massive. 
The massless case will be covered in the next session by taking the massless limit from the 
massive case. 

Now we illustrate the symmetric antler decay defined in Eq. ([7]), which has two indepen- 
dent rapidity parameters r/s and rj a : 



cosh r]B 



m D 



2m 



cosh rj a 



in 



B 



m 2 x + m 2 a 



B 



2m a m B 

BiB 2 decay, and rj a from B\ 



m n 



Note that r\B is determined from D 
B 2 — > CL2X2 equivalently). 

In the D rest frame, the momenta of the particles a\ and ai are 

^ cosh r]B cosh r\ a + sinh rjB sinh r\ a cos Q\ \ 
sinh r] a sin 9\ 


\ sinh rjB cosh r\ a + cosh r/s sinh r\ a cos Q\ ) 

( cosh r]B cosh r] a + sinh r/B sinh r] a cos 62 \ 
sinh rj a sin 9 2 cos 
sinh rj a sin 9 2 sin 
y — sinh r/B cosh r\ a — cosh r/s sinh r\ a cos 9 2 ) 

where the internal phase space angles of 9 2 , (p) are defined in Fig. |3J 



(17) 

aiXi decay (or 



(18) 



k 



m n 



(19) 
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(a) (b) 

FIG. 4: (a) The invariant mass m as a function of cos#i and cos#2> with = 0. (b) The 
normalized differential decay rate as a function of m. We take the mass parameters oim£> = 1 TeV, 
m B =400 GeV, m x = 200 GeV and m a = m Z - 



A. Invariant mass distribution 



The invariant mass of the two visible particles a\ and a 2 is explicitly obtained from k l 
and k ( 2 D) in Eqs. (fl8|) and (fT9l): 

{2 coshes cosh?7 a + sinh?7 B sinh 7/ a (cos 6>x + cos 2 )} 2 



(D) 



m 



— sinh 2 r/ a (sin Q\ + sin #2 cos 4>) 2 — sinh 2 r; a sin 2 9 2 sin 2 

— cosh 2 sinh 2 7/ a (cos 9\ — cos 2 ) 



(20) 



In Fig. H](a), we show the invariant mass m as a function of cos 9\ and cos 6*2- For the sake 
of illustration, we take = 1 TeV, m# = 400 GeV, m a = m^, mx = 200 GeV, and fixed 
= 0. The mapping of this non-trivial hypersurface onto the m yields a singular structure 
in the dT /dm distribution as in Fig. ID^b). To understand how this distinctive feature occurs, 
we study this mapping by examining the following some critical points: 

• Point (i): (cos#i,cos0 2 ) = (1, 1) 

Since a\ and a 2 move back-to-back in the D rest frame, their invariant mass becomes 
maximum. The rapidity of a% in the rest frame of D is the same as that of a 2 , such 
that I rji^ I = 1 77^ 
is 77^ l} = 2(t] B + rja). 



— t]b + rj a . Therefore, the relative rapidity of a 2 with respect to a x 
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• Point (ii): (cos0i,cos0 2 ) = (±1,=f1) 

One visible particle, say a±, moves in the same direction of its parent B\ with the 
rapidity of rji^ = T)b + Va, an d the other visible particle a 2 moves in the opposite 
direction of its parent with rji^ = \vb~ Va\- If Va > Vb, the directions of a± and a 2 in the 
D rest frame are the same, which implies V^ = T$-r{£\ If Va < Vb, the direction of 
ai and a 2 are opposite so that ^i" = r/if^ +rj^\ Therefore, regardless of the ordering 
of 7] a and 7] B , rjat = 2r] B - Note that two configurations of (cos^x, cos^) = (1,-1) 
and (cos #1, cos 6 2 ) = (—1, 1) are symmetric to each other. 

• Point (hi): (cos6>i, cos^) = (-1, -1) 

ai and a 2 move in the opposite direction to Bi and B 2 in their parent's rest frames, 
respectively. Their rapidities are \y]a^\ = \r] a ^\ = |?7b— ?7 a |, leading to rj^ 1 ' = 2\t] B — r) a \. 

• Point (iv): 8 2 = #i, = 0, cos^i = — tanh?7 B / tanhr/ a with r] a > t] B 

This special configuration gives rise to the same four-momenta of the two visible 
particles as can be seen in Eqs. (fl8l) and ( fl9l) . a\ and a 2 are relatively at rest, resulting 
in = 0. The condition 7] a > rjB is required to guarantee the equality of k[ D ^ and 
k 2 D \ which cannot be achieved if the particle B is boosted more highly than the 
particle a (or equivalently | cos^l < 1 for physical configurations). 

Point (i) corresponds to the maximum endpoint, and Points (ii) to the cusped peak. 
When Point (iv) exists, Point (hi) corresponds to the non-smooth kink, and Point (iv) 
to the minimum endpoint at m — 2m a . Otherwise, Point (hi) becomes the minimum 
endpoint at m = 2m a cosh(?7 B — r] a ). 

Now we present the analytic expression of the invariant mass distribution. The functional 
forms are different in the following three mass regions: 



In Fig. [51 we show the invariant mass distribution dT / d(m j Vn max ) for TZi, TZ 2 , and 
7?-3. Regardless of the parameter regions, the maximum endpoint of the m distribution 
corresponds to POINT (i): 



Ki ■ Vb < — 




< Vb < V 



Tli-Va< Vb- 



(21) 



m 



max 



2m a cosh(77 B + r) a ). 



(22) 
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d(m/m max ) 
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FIG. 5: The shapes of the function of dT/dm for three representative regions, 1Z±, IZ2, and IZ3. 
With the fixed r\ a = 1, we set t\b = 0.3 for TZi, rjs = 0.7 for IZ2, and tjb = 1-5 for 7^3. 

For 72 1 and 72.2, the minimum endpoint occurs at m = 2m a while for 72-3 the minimum is 
different: 



The condition of tjb > Va m 72-3 does not allow the equality of k\ = k 2 which would 
lead to m min = 2m a . In TZi and TZ 2 , there are two non-smooth points in the middle of the 
distribution. Let us call the point at the smaller value of m (marked by squares) a knee 
point and the other point at the larger value of m (marked by circles) a cusp point. In IZi, 
the knee point corresponds to Point (11) and the cusp point to Point (in). In 72-2, it is 
opposite. In 72-3, there is only one sharp peak, the cusp. We summarize the results of the 
minimum, cusp, knee, and maximum of the m distribution in Table [H 
The invariant mass distributions for three mass regions are 




2- 



(23) 




otherwise , 






(24) 
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K x 


K 2 


^3 


^min 


2m a 


2m a 


2m a c rjB - rja 




2Tfl a C^ B 


2m a c rjB ~ ria 




^cusp 


2m a c rjB ^ rja 


2vn a c r ^ B 


2ui a c r j B 


^max 


2m a c r]B+Va 


2m a c rjB + ria 


2m a c r)B + r)a 



TABLE I: Summary of the minimum, cusp, knee, and maximum of the m distribution for the mass 
parameter regions TZ\, IZ2, and R3. We have used a concise notation of coshx = c x . 



Region 


m D 


m B 


m a 


m x 


VB 


Va 


Mass-1 


Tlx 


650 


300 


m z 


100 


0.41 


1.06 


Mass-2 




850 


330 


m z 


100 


0.74 


1.18 


Mass-3 


n 3 


1000 


250 


m z 


100 


1.32 


0.80 



TABLE II: Test mass spectrum sets for the symmetric antler decay. All masses are in units of GeV 
and mz is the Z boson mass. t\b and r\ a are the rapidities of the particle B and a in its parent 
rest frame, respectively. 



dV 

dm 



2m cosh- 1 



oc < 



K 2 



m 

m 
0. 



2(r/ a -r/ B ) + cosh- 1 
2(r/ a + rj B ) - cosh" 1 



if 2m a < m2m a c VB _ Va ; 
if 2m a c VB - ria < m2m a c VB ; 

if 2m a c VB <m< 2m a c VB+ria , 
otherwise , 



(25) 



dr 

dm 



oc < 



K 3 



m 

m 
0. 



2r\ a -2t] B + cosh 1 



2 

m 

2m 2 „ 



2 Va + 2 VB - cosbT 1 



if 2m a c VB - Va < m < 2m a c VB ; 

if 2m a c VB <m< 2m a c 7lB+Va ; ( 26 ) 
otherwise . 



Here we have employed the narrow width approximation and ignored spin correlation effects. 
The detailed derivation is presented in Appendix 

In order to show the characteristics of the m distribution, we take three samples for mass 
parameters in Table [TT1 We label them as Mass-1, Mass-2 and Mass-3, each of which 
corresponds to the kinematical regions of TZi, IZ2 and TZs, respectively. The visible particle 
is assumed to be the Z boson. 
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FIG. 6: The normalized invariant mass distribution dT/dm for test mass sets in Table ITT1 

In Fig. [61 we present the invariant mass distributions for the mass parameters in Table 
ITT1 All three mass sets yield sharp cusp structures. The m minimum for TZi and 7^.2 is 2mz 
as discussed before. For the IZ3 case, however, fast-moving intermediate particle B yields 
^min = 2mzcosh(?7B — Va)- Unfortunately, we still have a two-fold ambiguity between IZi 
and 7^-2 because we do not know a priori whether the observed m CUS p is 2m a coshes — Va) or 
2m a cosh t)b- As shall be shown in the next section, the transverse momentum distribution 
breaks this ambiguity through its cusp and endpoint structures. In addition, the 1Z± and 
7^-2 cases have the knee structure, even though it is challenging to probe with the expected 
statistics at the LHC. 
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1 

Mass- 
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0.007 




/ \ Antler m a = m,z - 


0.006 
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0.004 
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0.002 




/ \ X. >^ Mass-3 - 


0.001 
n 


1 
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B. Transverse momentum variables: itit, Pt, and pxi 

In this section, we investigate the distribution of the transverse mass mj>, the transverse 
momentum variables px and pxi- In Fig. we show the my distribution. All the raj- 
distributions for TZ\, IZ2, and IZ3 do not have any cusped peak. The maximum in the nij 
distribution is the same as the maximum of m: 

( m T) max = m max . (27) 

The confirmation of the same maxima in the m and distributions will help the recon- 
struction of the antler decay. 
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FIG. 7: The normalized transverse mass distribution dT/drnx for test mass sets in Table ITTl 
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FIG. 8: The normalized transverse momentum distribution dT/dpT and dT/dpTi for test mass sets 
in Table M 

In Fig. El we plot the distribution of pr and pn- The total px distribution does not reveal 
the cusp structure, as expected from the rriT distribution. In addition, its maximum is at 
the end of a long tail, which is statistically disadvantageous to observe. The cluster transfer 
mass rric of ax-a^-fr system has no cusp structure either. 

The transverse momentum of "one" visible particle shows quite different distribution. 
First, note that one unambiguous px% distribution can be constructed out of two visible 
particles, because of the symmetric topology of the antler decay. This px% distribution 
shows the cusp structure as well as the fast-dropping maximum structure. The cusp and 
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FIG. 9: The normalized cos distribution for the massive visible particle cases, 
maximum of p?i are 

(PTi)cus P = m a \sinh(r) a -r) B )\, (28) 
(Pri) max = m a sinh(?7 a + 77 S ). 

Note that (pri) max gives the information about r] B + r) a , which is the same from m max 
in Eq. (1221 . Remarkable is that (pTi) cusp is common for all three regions 7^2,3, which 
determines \t]b — r) a \. By comparing (pri)cusp with m cusp , we can distinguish 1Z\ from 1Z2- 
This breaks the two-fold ambiguity in the measurement of m cusp for TL\ and "K.^- 

C. Angular variable: cosG 

We consider the distribution of cos0 defined in Eq. (I15jl . Here is the angle of one 
visible particle with respect to the cm. moving direction. As in the pxi distribution, the 
symmetric decay chains of the antler decay guarantee one unique cos distribution as shown 
in Fig. [9j All cos© distributions for TZi, Tl-2 and 7?-3 are symmetric about cosO = 0, and 
have sharp cusps. 

IV. MASSLESS VISIBLE PARTICLE CASE 

Now we consider the massless visible particle case. As suggested in Eqs. (J2J) , (SI) , and ([5l), 
many new physics processes for the antler decay have massless visible particles. Although 
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we cannot directly apply the results with the massive visible particle to this case since the 
rapidity t] a diverges, we can obtain the massless limit by using the finite combinations of 



m a c Va and m a s Va : 



.. ,. m B mi- \ , . 

hm m a c Va = hm m a s Va = —- 1 T . (29) 

A. Invariant mass distribution 

In the massive visible particle case, the functional form of the invariant mass distribution 
is different according to three mass regions of 1Z\, 1Z2-, and IZ3. In the massless visible 
particle case, only 1Z\ applies since r} B <S r\ a . Two locations of m m i n and mknee merge 
because m a = 0. The cusp and endpoints are given by 



m 



(0) 



0, (30) 
™ ( X = m B (l-^y^, (31) 

mSL = m B (l-^)e^. (32) 

V m B J 

Here the superscript (0) is used for emphasizing m a = 0. The product of the cusp and the 
maximum is 

™£WSL = < ( 1 - 4V , (33) 

V 171 B / 

which depends only on the second step decay of B — > aX. The ratio is 

(o) 

^W 2 -, (34) 

which is determined only by the first step decay of D — > BB. 
The invariant mass distribution is simplified into 



i {0) ■ 



dT 

dm 



m log (^^ff 1 ) ' if < m < 

\ ^cusp / 

m\og(^^j, if mffip < m < m£L; (35) 

0, otherwise. 



For < m < mcusp, dT /dm is a linear function of m. For mi% <m< m£L, it is a concave 
function with the maximum at m = mSL/e. Depending on the relative position of miulp 
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m D 


m B 


m x 


m a 


(0) 
^cusp 


(0) 
^max 


1 cos6| max 


Mass-1 


1000 


470 


440 





40.7 


82.9 


0.34 


Mass-2 


1000 


440 


410 





34.6 


97.1 


0.47 


Mass-3 


1000 


400 


370 





28.9 


115.5 


0.60 



TABLE III: Test mass spectrum sets for the symmetric antler decay with massless SM particles. 
All masses are in units of GeV. 
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FIG. 10: The normalized invariant mass distribution dT/dm in the massless visible particle case 
for the mass parameter sets in Table IIII1 

and mmL/e, the maximum of the concave function may or may not show in the function 
of dT/dm, which determines the sharpness of the cusp. If mSL/e < mSlp (or equivalently 
me > 0.443m/)), dT/dm is linearly increasing up to m = mSlp, and decreasing after that: 
the cusp is sharp. If mcul p < mmL/e, dT/dm keeps increasing after m = m™ sp , reaches the 
maximum of the concave function, and finally falls down: the cusp is not sharp. The m cusp 
structure is most useful when the D — > BB decay is near the threshold. 

In order to show the functional behaviors specifically, we take three mass sets for the 
massless visible particle case in Table HTT1 The mass parameters in the Mass-Io correspond 
to the case where both the first decay D —> BB and the second decay B — > aX occur near 
the threshold. This is motivated by the decay of the second KK mode of Z boson in the 
UED model in Eq.(j4]). The MASS-2o represents the marginal case for the sharp cusp, i.e., 
ms~ 0.44 mj). The Mass-3q case has large mass gaps. 
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FIG. 11: The normalized transverse mass distribution dT/drriT for the massless visible particles. 
The mass spectrum sets are described in Table ITO1 

Figure (TU] shows the m distributions. All three mass sets in Table II I II have the same 
^mm = 0. The sharpness of the cusp structure is different. The nearly degenerate mass case 
(Mass-1o) has a very sharp cusp. The marginal case (Mass-2o) shows also an observably 
sharp cusp. The large mass gap case (Mass-3o) has a rather smooth cusp. If the number of 
events is not enough, the obtuse cusp in the Mass-3o is difficult to read. The measurement 
of the m cusp is still possible since the functional form in Eq. ( |35l) can be used to fit the 
data and to read the cusp position. 

B. Transverse momentum variables: my, px and pxi 

Now we turn to the kinematic variables involving transverse monentum. First, the my 
distribution in the massless visible particle case does not show any cusp structure as shown 
in Fig. [TTJ The absence of my cusp is a common feature of the antler decay. The 
maximum stands at the end of fast-dropping function for all three mass sets, which is easier 
to read. In addition it is the same as the m maximum: 



Figure d2] shows the distributions of the total transverse momentum and individual 
PTi- As in the massive visible particle case, the total pt distribution is very smooth and 
gentle, without any cusp structure or fast dropping maximum. Instead, the pt% distribution 



= m. 



max' 



(0) 



(36) 
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FIG. 12: The normalized transverse momentum distribution dT/dpT and dT/dpTi for the massless 
visible particles. Test mass sets are in Table [TTT1 



shows very sharp cusp, much sharper in general than the invariant mass distribution. Even 
the Mass-3o case, which suffers from the dull cusp in the m distribution, has a very sharp 
PTi cusp. In addition the pt% maximum is at the end of a faster dropping function. 



The analytic expressions of (ptJcSL an d (P^JmL can ^ e easily obtained from Eq. ( 128]) by 



'cusp 



(0) 
max 



applying Eq. f f29|) : 



ffsp and (PTi) { 2 

which is another way to check the antler decay topology. 



'/max 2 max ' 



,(0) 



The measurements of (pri)^ sp and (pri)mL P rovi de the same information as mSsp and mmix, 



(37) 
,(°) 



C. Angular variable: cos© 

Figure [TBI shows the normalized dT/dcosQ distributions for three massless visible particle 
cases. The function increases with |cos0|, and drops to zero suddenly at |cosG|^ x . This 
is because the cusp and the endpoint merge, resulting in more pronounced endpoints with 
sharp peaks at both ends. The maximum of cos is simply determined by the first step 
decay D ->• BB: 

l cos0 lSL = tanhr?B- (38) 
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FIG. 13: The normalized cos distribution for the massless visible particle cases. 
The full analytic function of d$4/<icosO is given by 



The suddenly ending behavior of the cos B distribution is because massless visible particles 
cannot access all kinematic space of cos 9. The detailed derivation of Eqs. fl38j) and fl39|) is 
in Appendix IA 41 

V. EFFECTS FROM REALISTIC CONSIDERATIONS 

In the previous sections, we have considered the kinematics only, ignoring the decay width 
of the intermediate particle B, the longitudinal boost of the parent particle D, and the spin 
correlation. These S'-matrix element effects can smear the kinematic cusps and endpoints. 
In the following, we discuss the limitation of determining the missing particle mass using 
kinematic singularities. 

A. Finite width effects 

The previous results are based on the narrow width approximation. This approach is 
very effective for the proposed processes in the MSSM, Z' super symmetry, UED, and LHT 
models since all the intermediate particles (x^ ^ S > and t_) have very small total decay 
widths, much smaller than one percent of their masses. If the total decay width T# is 




otherwise. 




(39) 
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FIG. 14: The invariant mass distribution dT/dm in the antler decay with finite Tb effects. We 
have taken the mass spectrum sets of Mass-1 and Mass— lo- 

large, its effects can smear the cusp and endpoint structures. If the on-shell B particle 
is kinematically not accessible so that the decay process is through off-shell B, then the 
singular structures are destroyed completely since there is no constraints on the phase space 
from the mass relations. 

In Fig. [HI we show the invariant mass distributions with the effect of finite Tb for the 
massive SM particle case (Mass-1) and the massless case (MASS— 1 ). We take T B to be 
3%, 10%, and 50% of tub for the massive case, and 1%, 10%, and 50% for the massless case. 
If Tb/Mb is small enough (< 3% for the massive case and < 1% for the massless case), the 
m cusp remains fairly preserved. Even though the sharp cusp gets dull slightly, the position 
of the cusp is not shifted for both cases. The endpoint position is stable for the massive case, 
but shifted considerably for the massless case. If Tb/ttib is about 10%, the cusp is smeared 
into a round peak and the endpoint position is shifted significantly for both cases. Still the 
peak of the smeared cusp stands at the same cusp position. If Tb/Mb — 50% in which case 
a large contribution to the S'-matrix element arises from the intermediate off-shell B, the 
sharpness and position of the cusp are lost. The endpoints move towards new positions of 
m = tud — 2m>x- This is from the allowed phase space of the decay D — > XXaa. At least 
we can determine the mass difference between D and X using the m distribution. 

Now we show the Tb effects on the p?i distributions in Fig. [15j We take the massive 
Mass-1 case and the massless Mass-Iq case. The pn distribution has very vulnerable 
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FIG. 15: The normalized transverse momentum distribution dT/dpT with various finite total decay 
width of the intermediate particle B. 
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FIG. 16: The normalized cos distribution with various finite total decay width of the intermediate 
particle B. 



cusp and endpoint from the finite Tb effects. Even for small width effects (< 3% for the 
massive case and < 1% for the massless case) the sharp cusp becomes dull, and its position is 
significantly shifted. The pti maximum is more sensitive to the V g effects. For the massless 
SM particle case, even 1% ofTs/m^ shifts the position of p™f x a lot. 

Finally, Fig.[16]shows the cos distribution with finite Tb effects. Here the most dramatic 
collapse occurs. Even with very small width of T B /m B = 1%, the sharp cusp becomes round, 
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FIG. 17: Comparison of the cosO distribution in the D rest frame (thin curves) and in the pp 
lab frame with y/s = 14 TeV (thick curves). The mass parameters used here are Mass-1o and 
MASS-2 defined in Tables M and HD 

difficult to read. For Y B /m B = 10%, the cusp shape is lost completely. 

In summary, the effects of the finite width of the intermediate particle B smear the cusp 
shape and shift the cusp position to some extent. The invariant mass distribution has the 
least distortion, while the pt% and cos distributions have significant changes, especially for 
the massless visible particle case. However, the proposed processes in Eqs. (JT|) — (J4|) are not 
affected since Tb/tub is much smaller than 1%. 




B. Longitudinal boost effect 

At hadron colliders, the longitudinal motion of the particle D is not determined even when 
D is singly produced. Among the discussed kinematic variable, only the cos0 is affected, 
which is defined with a momentum in the D rest frame. In order to see the longitudinal 
boost effects, we convert the cos distribution in the D rest frame into that in the pp frame 
at the LHC, by convoluting with the parton distribution functions of a proton. We have 



used CTEQ6 25]. In FigJTTl we compare the normalized cos© distribution in the D-rest 
frame (thin curves) with that in the pp lab frame with y/s = 14 TeV (thick curves). For 
simplicity we assume that the heavy particle D is produced through the s-channel gluon 
fusion and/or qq annihilation. 
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FIG. 18: The invariant mass distribution dT/dm and the individual transverse momentum dis- 
tribution dT/dpTi with and without spin correlations of -> L^L^ -> i'B^h+B^ i n the 
minimal UED model. We have set 1/R = 500 GeV and AR = 20. 

In the massive visible particle case (mass-2), the cusped peaks vanish almost completely. 
In the massless case (mass-1o), the pointed cusps become round, very hard to read. We 
conclude that the cusp in the cos distribution is not observeble at the LHC. In the e + e~ 
collisions, however, the fixed cm. energy removes the longitudinal boost ambiguity, and 
thus the cos G cusp provides valuable information on the missing particle mass. 

C. Spin-correlation effect 

The effects of the spin-correlation by the full matrix elements are different from new 
physics process to process. In addition, if we consider the associated production of the par- 
ticle D in order to control the SM background, the spin correlation effects get intertwined 
with the additional px and/or longitudinal boost effects. To maximize the discovery signifi- 
cance, it is desirable to develop an individual strategy for each process in Eqs.([T])-(jl|), which 
is beyond the scope of this paper. 

Generically, the positions of cusps and endpoints are not affected by the spin correlation 
effects since they are determined purely by the constrained phase space, i.e., by the mass 
relations In order to see this feature, we consider the decay in the framework of 



26 



the minimal UED model (mUED) 261 ]: 



Z (2) L (D + L (D £- ^(1) + (40) 

In Fig. [TU we show the m and pTi distributions including the full matrix elements of 
the process in Eq. @Q} at the LHC with y/s = 14TeV. We have fixed 1/R = 500 GeV and 
AR = 20, which generates the KK masses of = 1048 GeV, rris = 515 GeV, m a = 0, and 
mx = 500.9 GeV. First finite width effects are negligible: very degenerate mass spectrum 
in the mUED model yields very small total decay width such that Tb/ttlb ~ 10 -4 . Second 
the longitudinal boost effects do not apply to m and pxi- As shown in Fig. [181 the spin 
correlations hardly change the m and pt% distributions. The distributions with and without 
the spin correlation effects are almost identical. 

Brief comments on the SM background and detector simulation effects are in order here. 



In Ref. [18], we have shown that the cusp structure survives over the SM backgrounds and 
the detector simulations in a benchmark process of pp — >■ Z' — > £ + £~ — > £ + Xi^~Xi in a 
supersymmetry model with an extra U(l) gauge field. In addition the missing particle mass 
as well as the intermediate particle mass can be determined, even though the uncertainty 
is about 10%. It was demonstrated that the analytic expression for the invariant mass 
distribution is very helpful to reconstruct the mass parameters by best-fitting. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we have considered the antler decay topology of a parity-even heavy particle 
into two missing particles (X\ and X 2 ) and two visible particles (a\ and a 2 ) via intermediate 
on-shell particles (Bi and B 2 ). We studied the singularity structures in various kinematic 
distributions, especially non-smooth peaks called the cusps. We show that the distributions 
of the invariant mass m of ai and the individual transverse momentum pxi, and the 
cos develop conspicuous cusp structures. We have provided the detailed derivations for 
the positions of the cusps as well as the endpoints in terms of the particle masses. The 
analytic functional forms of the invariant mass and cos distributions have been also given. 

The cusp and endpoint structures of the antler decay have a few advantages: (i) if the 
parent particle mass is known from other decay channels, they can be used to determine 
both the missing particle and intermediate particle masses; (ii) the cusped peaks are more 
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identifiable than endpoints and kinks due to higher statistics at the kinematical maxima; 
(Hi) the simple configuration of outgoing particles, two visible particles and two missing 
particles, avoids combinatoric complication, which is troublesome in many missing particle 
mass measurement methods; (iv) the position of the cusp is independent of the S'-matrix 
element such as the spin correlation effects, since it is purely determined by the phase space. 

We point out that the pTi cusp and endpoint have some desirable features for observation. 
The pTi cusp tends to be sharp irrespective of mass parameter regions. It is complementary 
to the robust m cusp, which is sharp only when the masses are nearly degenerate. The pTi 
endpoint is always located at fast-dropping end, which is easier to read off. Finally, the cusp 
position for the massive visible particle case is uniquely determined by the involved masses, 
while the m cusp has two-fold ambiguity. 

It is noted that the cusp structures have some limitations for the missing particle mass 
determination, especially at the LHC. The sharp cusped peaks in the cos distribution are 
not readily observable at the LHC, due to the longitudinal boost of the produced D particle. 
The effects of the finite width of the intermediate particle could affect the cusp and end- 
point in the individual transverse momentum distribution. However, for generically weakly 
coupled theories beyond the SM, the new particles for the antler decay have relatively small 
decay widths, and thus the pxi cusp is expected to be preserved. The cusp in the invariant 
mass distribution, which is the most robust observable at the LHC, is most pronounced for 
a degenerate mass spectrum. 

In addition, the relations among different cusps and endpoints help to identify the antler 
decay topology. For example, the my maximum is equal to the m maximum. The cusp 
and endpoint of pt% distribution are half of those of m distribution in the massless visible 
particle case. One can use these facts for the consistency of the assumptions on the event 
topology. Similar intriguing relations exist for the massive visible particle case. 

In conclusion, if a new physics model accommodates an antler decay, the measurement of 
kinematic cusps and endpoints can be helpful to determine the missing particle mass as well 
as the intermediate particle mass. The proposed processes in various new physics models are 
expected to have stable cusp and endpoint structures in the m, my, and pt% distributions 
at the LHC. 
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Appendix A: Kinematic distributions in the symmetric Antler decay 

In this Appendix, we derive the kinematic distributions of the invariant mass and the 
cos for the symmetric antler decays D — > B\ + B 2 — > a\Xi + a 2 X 2 where =m j e 2 =ms, 
m ai = m a2 = m a and mx 1 = f«i 2 = rrix- We first describe the four-body phase space and 
choose special internal phase variables. For each four-momentum, we specify the reference 
frame in the superscript. For example, k[ Bl ' is the four-momentum of ax in the rest frame 
of B\. Since we calculate physical quantities mostly in the D rest frame, we omit the 
superscript for the D rest frame for simplicity. Note that this is different from the notation 
in the main text where we omit the superscript for momenta in the lab frame. 



1. Four-body phase space 



We consider four body decays of 



D{P) -> ai(fci) + a 2 {k 2 ) + Xi(fc 3 ) + X 2 (k 4 ). 



(Al) 



The differential decay rate of the process is 

1 



dT 



2m 



D 



\M\ 2 d<S> A (P; k x ,k 2 ,h,h) 



(A2) 



where \M.\ 2 is the he 



space, defined by 



27 



icity amplitude squared, and <i$4 is the element of four-body phase 
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(A3) 



If the decay in Eq. (IA1I) is through the antler decay, i.e., through D — > B\B 2 followed by 



Bi — > diXi(i = 1,2), the helicity amplitude squared \A4\ 2 has two propagator factors of B\ 
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and .E?2- Using the narrow width approximation Tb/tub <C 1, the matrix element squared 
can be expressed in terms of two Dirac delta functions: 

1 1 



l " M ' 2 ~ lM? (Pi - m%y + m%T% b|-m|) 2 + m|r| 

r S «- B ^ 2 (_^_y ^ - m|) 6(pl - m B ), (A4) 

where pi = k\ + and p 2 = ^2 + ^4 are the momentum of B\ and _E>2, respectively. In this 
limit, \M.\ 2 does not develop any singular behavior and still remains as a smooth function 
containing spin correlation information. 

After the integration using delta functions, the differential decay width is simplified to 

1 — ■ y /2 

dT = — — — ^—75- \M. I 2 A R \ a d cos 9id cos 9 2 dd), (A5) 

2 15 ir 4 m D m B T B 

where A# = A (1, m B /m 2 D , m 2 B /m 2 D ), X a = X(l,m 2 /m 2 B ,m 2 K /m 2 B ), and the standard kine- 
matic function is A(a, b, c) = a 2 + b 2 + c 2 — lab — 2ac + 2bc. The polar angles of Q\ and #2 
and the azimuthal angle <fi are defined in Fig. |3j For simplicity, we use short-hand notations 
of 

V\ = cos#i, u 2 = cos# 2 > (A6) 
and name dv\dv2d(f) the normalized four-body phase space <i$4 of the antler decay: 

d$4 = dv\dv2d(j). (A7) 

2. Change of variables and the independence between angular variables 

For a general two body decay of a — > be: 

a(Pa) -> KPb) + c(p c ), (A8) 
the energy-momentum conservation in the rest frame of the parent particle a leads to 



eW+E®, (A9) 



p£ o) = -Pi a) - 
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From the on-shell conditions of p 2 = m 2 (i = b,c), the energies and momenta of the particles 
b and c are simply expressed by the rapidities rjb and rj c : 



E 



E. 



(a) = m 2 a + m 2 b - m 2 c 
2m a 

(a) _ m 2 a -m 2 + m 2 c 



2m n 



mb cosh rjb, 
m c cosh rj c , 



Pi 



(a) 



mb sinh rjb = m c sinh rj c 



(A10) 

(All) 
(A12) 



For the symmetric antler decay, the same masses of ms = m Bx = mB 2 and m a = m ai = m a2 
lead to two independent rapidities: 



mo rn 2 R — m\ + m 2 

cosh?7 B = , cosh?7 a = 



2m 



B 



2m a m B 



(A13) 



Now we present less intuitive but more convenient kinematic variables. First, we consider 

( B ) 

the rapidity of 02 in the rest frame of B\, not B2, denoted by a = rja 2 . 



cosh a = cosh2?7B cosh?7 a — v 2 sinh.2rjB sinh?7 a , 



(A14) 



where v-i is defined in Eq. ( 1A6I) . The second useful variable is u, the cosine of the angle 6 ai 



(Bi) 



(12 



between a\ and a 2 in the rest frame of B\. 



u 



k (Bi) . k (Bi) 



(Bi) 



SBi) 



I — vfy/l — v\ cos0 + (sinh2?7 B — v 2 cosh.2rj B ) v 1 
cosh2r/ B — v 2 sinh2?7 B 



(A15) 



For simplicity, we define 



cosh2?7 B — v 2 sinh2?7 S , 
sinh2?7 B — v 2 cosh2rj B . 



(A16) 



Then the azimuthal angle <fi is inversely obtained by 



COS ( 



uv' 2 — V1V2 



(A17) 



The advantage of this new angular variable u is that d 2 ^^/ dudv 2 = n: u and v 2 are 
independent variables contrary to the expectation from the functional dependence of u on V\ 
and v 2 in Eq. (IA15|) . In order to show this non-trivial result, we begin with gP$ 4 / 'dvidv 2 d(p = 
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1 in Eq. (1A7j) . We change the variable <fi into u as 
o?$4 = dv\dv 2 d<\> = dv\dv 2 du 
= dv\dv 2 du 



du 



(A18) 



v' 



dv\dv 2 du- 



v' 



y/f(u,V 1 ,V 2 ) 

Since the integrand v' 2 j v±, v 2 ) is not separable into products, u, v\ and v 2 are not 
independent with one another. If we integrate one of the three variables, however, we have 
the statistical independence of the remaining two variables. First v i and v 2 are independent 
variables by definition. In order to see the independence of v 2 and u, we integrate v\ out for 
given u and v 2 . The integration limit of v\ is matched with the roots of f(u,vi,v 2 ) = for 
fixed u and v 2 . The result of the integration is a simple constant: 

(max) 

dvi , V2 = 7T. (A19) 



" nin> y^i-^) (^! max) -- 2 ) w 



Therefore d^^/dudv 2 = n is also flat: u and v 2 are independent. Similarly, one can show 
the independence of u and v\ from the symmetry under the exchange of v\ and v 2 . 

3. The invariant mass distribution 

The invariant mass m of a\ and a 2 is more simply expressed in terms of a and u by 

m 2 = 2m 2 a + 2m 2 (cosh ?7 a cosh a — u sinh ?7 a sinh a), (A20) 

where a and u are defined in Eqs. (1A14I) and ( 1A15I) respectively. The expression in the 
parenthesis of Eq. ( 1A20I) is nothing but the cosine hyperbolic of the rapidity of the particle 
a>\ in the rest frame of a 2 : 

t 
2m 2 a 

Now let us change variables from (u,v 2 ) to (x,®)'- 

dudv 2 = d\da k — . (A22) 

sinh 2t]b sinh t] a 



cosh = - — - — 1 = cosh rj a cosh a — u sinh r] a sinh a. (A21) 

a 



32 



cosh(2?] 1 + 2-q B 



cosh(2r/i — 2t]b 



cosh2r/£ 




la 

a 



Va + 2ve 



FIG. 19: Allowed parameter space of (a,x) plane for the region TZi. 
Note that the Jacobian factor is simply a constant. From <i 2 $4/ ' dudvi = vr, we have 

d$A IT /-<W(X) 



da, 



(A23) 



dx sinh2?7 B sinh r\ a j arnin (x) 
where a m i n (x) and a max (x) are the minimum and maximum of a variable at a given x, 
respectively. 

In order to obtain Q: m in(x) and a max (x)) we use the conditions of u G [—1, 1] and t>2 G 
[—1, 1]. Then the definitions of cosh a and x m Eqs. ( 1A14I) and (1A21I) . respectively, constrain 
the values of cosh a and \ as 



cosh(2?7 B — r] a ) < cosh a < cosh(2?7 B + r] a ) 
cosh(a — T) a ) < x < cosh(a + r) a ). 



(A24) 
(A25) 



Therefore, the values of a m in(x) an d «max(x) i n Eq. ( 1A23I) depend on the relative size between 
rjB and rj a /2 or rjs and r) a . This is related with the three different mass parameter regions 
of Hi, 7Z 2 , and 7Z 3 in Sec. IIII A[ 



(A26) 



i<-i ■ Vb < y, aC 2 : — < Vb < Va, TZ 3 : r] a < 7} B . 

Let us elaborate the derivation of a m in(x) an d «max(x) f° r the region H 1 . Figure 
illustrates two curves of x — cosh(a — r) a ) and x — cosh(a + r) a ) in the parameter space of 
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(a, x)- Within the bound of rj a — 2rjB < a < r] a + 2t]b as in Eq. (1A24p . a m in(x) an d «max(x) 
are different according to the value of Xi summarized by 



For Tlx 


< -^min(x) 




J da 


1 < X < c 2vb 


rj a - cosh" 1 x 


rja + cosh" 1 x 


2 cosh" 1 y 

( 


C2 VB < X < C2r) a -2 VB 


Va ~ %Vb 


Va + 2?7b 




C2 Va -2 VB < X < C2r lB +2r la 


-rj a + cosh" 1 x 


Va + 2f] B 


2r] B + 2t] a - cosh" 1 x 



(A27) 



Here we use the simplified notation of c x = coshx. The derivations for 7?. 2 and TZ^ are 
similar and straightforward. 

With the help of Eq. (1A23[) . the final expressions for d&±/dm is given by 



1 d$A 



N dm 



2m cosh 1 ( 



if 1 < X < c 2vB ; 

^Vsm, if c 2vb <X< c 2( Va - VB y, ( A28 ) 

^ 2{j] a + T) B )m - m COsh" 1 (j^ - 1 J , if C 2 (r, a -r, B ) <X< C2( Va + VB ), 



1 d$* 



N dm 



n 2 



1 d&i 



N dm 



K 3 



2m cosh" 1 (jgg - l) , if 1< x < c 2(f?a -r 

2(?7a - VB)m + m cosh" 1 (j^ - 1 V if c 2 ( J7a _ ?7s ) < x < c 2f?B ; (A29) 

2(% + 7?b)™ - m cosh" 1 - l) , if c 2r , B <X< c 2 (r, a + VB ), 

2(Va ~ Vb)™ + m cosh" 1 - 1 
2(Va + Vb)™ - m cosh" 1 - 1 



if C2( Va - VB ) < X < C 2 VB ] 

If C 2r , B < X < C 2(Va+VB)l 



(A30) 



where the normalization factor A" is 

A 



TV 



1 



smh2r]B (m a sinh?7 a ) 2 
4. The angular distribution dT/dcosQ 



(A31) 



In this subsection, we derive dQ^/dcosO, restricting ourselves to the massless visible 
particle case (m ai = m a2 = 0). Recall that G is the angle of a visible particle, say a 1; in the 
cm. frame of a± and a 2 , with respect to their cm. moving direction in the D rest frame. 
For d^i/dcosQ, we begin with <i 3 $4/ ' dv\dv 2 d(\) = 1. The key point is the Jacobian factor 
from dv\dv2d(j) to g?cos0. For this goal, we first obtain the analytic expression of <fi in terms 
of vi, v 2 , and cos 6. 
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In the case of m a = 0, the k\ and k 2 four-momenta in the B\ rest frame become 



fc (*0 = E [ Bl \l,Jl-viO,-v 1 ) = E e (l,k 1 



(A32) 



fc (B0 = s (B l)(u / J 1 _ v 2 CQS(f) Jl- V *s in< f > ,- V Z)=Etv' 2 0.,k 2 ), 



where Eg = tob(1 — m 2 K /m 2 B )/2, kj = kj/|kj|(i = 1, 2), and the definitions of v 2 and v 2 are 
in Eq. (IA16|) . Defining k = k± + k 2 — (E cra , k), we have some useful expressions of 



m 2 = 2Ejv' 2 {l -u), 

= E e (l + v> 2 ), 
>| 2 = E 2 i {l + 2v' 2 u + (v' 2 ) 2 }. 



(A33) 



\k^ 



Now the Lorentz transformation matrix from the B\ rest frame to the cm. frame of a\a 2 



is 



\ 



\ 



7cm 
k(Bl) 

m 



m 



/ 3 x3 + (7c m -i)k^)k^) r J 



(A34) 



where 7 cm = Ecm 1 /m, and the superscript T denotes the transpose of the vector. The 
three-momentum of a\ particle in the cm. frame of a\ and a 2 is 

E {Bi) 



(cm) 



m 



-k^) + k 1 + (7 cm -i)(k^)-ki Bl) )k^) 



(A35) 



Since k in the D rest frame is k^ = -P£ m) = -A^ a2 ^ B ^P^ l] , we have 



E 



(Bt) 
D 



-k + P 



m 



(Si) 
D 



( 7cm - 1) ■ p(*')k(*) . 



(A36) 



The dot-product of k[ cm ^ and leads to cos0: 



cos 6 



k (cm) . k ( D ) 



k; cm) iik^)| 



(A37) 



Finally we express cosB in terms of (vi,v 2 , 



cos 6 



lv 2 -Vi)S 



'ID 



^2 - \{v± + i; 2 ) 2 + \{ Vl - v 2 ) 2 c 2vb + 2v^T^ ^)(1 - C os0 



(A38) 
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Note that the maximum of cosG occurs when v\ = ±1 and v 2 = =Fl, i.e., when the visible 
particles a\ and a 2 are moving in the same direction. The maximum of cos in the -D-rest 
frame is then 

l cos0 lma X = tanhr/ B . (A39) 
Finally cos0 is expressed in terms of v±, v 2 , and 0: 

cos </> = \ I , (A40) 

where s v = sinh?7 for simplicity. For the Jacobian factor, we introduce three independent 
variables, v+, u_, and t, defined by 

2s 2 

v± = V! ± v 2 , t = — f- - c 2m . (A41) 
cos^ 



Note that the maximum of | cos 0| in Eq. (IA39|) leads to the t integration range as 1 < t < 00. 
Since 



<i$4 = dv\dv 2 d(j) = -dv + dv-dt 
the differential four-body phase space with respect to t is 

4 1 / , j2 1 



Of 



(A42) 



- / dv+dv 2 , (A43) 

dt A J ~ - W) + St - 2v'it -8-2v 2 + V ! 

The integration range is 

o -8 - 2vl + 8t - 2vl t , K . 

< vl < ± —, (A44) 



t 2 - 1 



, / 4(t - 1) 



< v j < -1 L. (A45) 

~ + ~ t+1 K J 

Finally the integration over V- and v + yields 

= aV2tt sinh 2 Vb—^t. (A46) 
d cos sin 15 

Appendix B: The invariant mass distribution of generic antler decays 

In this section, we present the analytic expression of the invariant mass distribution of 
generic non-symmetric antler decays with 7^ rriB 2 -i m a 1 7^ m a 2 an d mx-i 7^ m x 2 - The 
derivation is very similar to Appendix |AJ but in this general case the mass parameter space 
is divided into finer twelve regions. Since the derivation of the formulae for each region is 
long and tedious, we show only the results here. 
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1. Massive visible particles (m a ^ 0) 



In generic antler decays, there are in general six different rapidity parameters, given by 



cosh 7]x 1 
cosh rjB 1 
coshr/ ai 



Bi 



m 



m 



<'i 



2m Xl m Bl 



m 



D 



m 



m 



2m Bl m D 



m 



Bi 



m 



«i 



m 



-Vi 



2m ai m Bl 



COsh 7]x 2 

cosh r]B 2 
cosh r] a2 



m B 2 + m X 2 - m l 2 



2m X2 m B2 



m D + m B 2 



m 



Bi 



2m B2 m,D 

m B 2 + m l 2 ~ m X 2 



2m a2 m B2 



We define 



V++ = VB! + Vb 2 + Vai + r]a 2 , 

V+- = \VBi + VB 2 + -Va 2 \, 

V-+ = \VBi + VB 2 ~ Vai +Va 2 \, 

V— = \VBi +VB 1 - Vai ~Va 2 \- 



(BI) 

(B2) 
(B3) 
(B4) 
(B5) 



From positive definite definition of the rapidity, v++ is the larges among four 77±±'s. However 
the relative size of the other three 77's is different according to the mass parameters. We 

order rj+-i V-+ an d V an d name them to be 771 < 772 < Vs- We have 6 regions depending 

on this ordering: 

V+- < V-+ < V— » V-+ < V+- < V—> 
V+- < V— < V-+, V-+ < V— < V+-, 
V— < V+- < V-+, V— < V-+ < V+- 

To obtain dT/dm, we introduce the general x, defined by 

22 2 
v m — mi — rr 

J ai 

la 2 



X = cosh 77^"^ 



2m ai m a2 



(B6) 



The general invariant mass distribution have 12 different cases in total, given by 

• ^ \VBi + VB 2 ~Va 2 \> Vai Or T) Bl + Tfe + 7] a2 < Tj ai , 



1 d$ 4 
N dm 



- Vl m + m cosh" ( 2m ^ ma2 ° 2 ) , if c m < X < c V2 , 



V2 - Vi, 

V+ 

0. 



cosh 



22 2 

_l 1 m —m„_, —mi 



2m ai m a2 



if Cjj 2 <J X — Cf]3i 
if c^g < X — C V++ 1 
otherwise. 



(B7) 
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/ / 2 2 

1 / m —m?,^ —m* 



If \VB 1 + VB 2 ~Va 2 \< V ai < VB 1 + VB 2 + Va 2 , 

2m cosh 

1 



N dm 



2m ai m a2 

—r]im + mcosh 
(?7i + i] 2 )m, 
rj ++ m — m cosh 

Here the normalization factor N is given by 



2m ai m a2 



2m ai m„ 2 



N 



TV 



if 1 < x < c Vl , 
if c m < x < c^, 

if C»y 2 ^ X — ^773, 
if c r; 3 ^ X — C V++ ' 

otherwise. 



(B8) 



(B9) 



m ai m a2 sinh 277 B sinh r] ai sinh 77 fl2 
Note that the minimum of the invariant mass distribution can be different from m ai + m a2 , 
according to the mass parameter regions. Crucial is whether the kinematic configuration 
that a\ and a 2 are relatively at rest is allowed. 



2. Massless visible particles (m a = 0) 



In this subsection, we present the invariant mass distribution for massless visible particle 



but different intermediate particle cases, i.e., when ^ m B2 and m ai 



m 



a 2 



0. In 



this case, t] is always larger than 77 + „ and rj (-, leading to 773 = 77 Here we need to 

consider only the leading terms of O (nT>ai m a 2 ) , which are absent in cosh?7 + _ and cosh?7_ + . 
Therefore, the invariant mass distribution is divided into three regions. Using cosh" 1 x = 
ln(x + v 'x 2 — 1) f» ln(2x) for x > 1, we have 



where 



d$4 
dm 



oc < 



m log 
m log 
0. 



(0) 

1 

■ 

(0) 



(0) 



if < m < m^ 



cusp 1 



■f (0) ^ ^ (0) 

it m cu kp <m < mw; 
otherwise, 



(0) 
cusp 




rn 



(0) 



m 



Bi 



m 



m Bl 



m 



Bo 



m 



x 2 



m B2 



exp 



exp 



_Vb 1 ±Vb 2 
2 

rjBr + Vb 2 



(B10) 



(Bll) 
(B12) 
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This is the generalized results of Eqs. ( I31~j) and fl32l) . Note that the product mc Usp ?Timax 
depends only on the second step decays of Bi — >■ X^ai and B 2 — > X 2 a 2 while the ratio 
^max/^cusp only on the first step decay of D — > BiB 2 . 
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